Random walks on complex trees 
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We study the properties of random walks on complex trees. We observe that the absence of loops 
reflects in physical observables showing large differences with respect to their looped counterparts. 
First, both the vertex discovery rate and the mean topological displacement from the origin present 
a considerable slowing down in the tree case. Second, the mean first passage time (MFPT) displays a 
logarithmic degree dependence, in contrast to the inverse degree shape exhibited in looped networks. 
This deviation can be ascribed to the dominance of source-target topological distance in trees. To 
show this, we study the distance dependence of a symmetrized MFPT and derive its logarithmic 
profile, obtaining good agreement with simulation results. These unique properties shed light on 
the recently reported anomalies observed in diffusive dynamical systems on trees. 

PACS numbers: 89.75.Hc, 05.40.Fb, 05.60.Cd 
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I. INTRODUCTION 

Diffusion problems on tree structures pop up in a 
wide range of scientific domains, such as theoretical 
physics [l[ 0, Q, computer science [3, phylogenetic 
analysis [H and cognitive science [3]. Moreover, dynam- 
ics in tree structures have gained a renewed interest in the 
physics community as a spin-off of the attention devoted 
to the structural properties of complex networks M, M 
and dynamical processes taking place on top of them [l0( . 
Thus, along with the widely explored scale-free (SF) net- 
works [ll| , also SF trees have started to be used as under- 
lying topologies for dynamical processes. Interestingly, 
the absence of loops in trees turns out to have a strong 
impact on the considered dynamics, and relevant differ- 
ences between looped networks and tree topologies have 
been recently reported in several dynamical models, such 
as the voter model [l2[, the naming game fl3l. the ran- 
dom walk and the pair-annihilation processes [14], and a 
model for norm spreading [ll|. The properties of most 
dynamical processes on looped network can be reason- 
ably accounted for by annealed mean-field theories [Tol ] . 
which rely only on information about the degree distribu- 
tion and degree correlations [lfj], and consider the net- 
work as maximally random at all other respects. The 
behavior observed in trees, different from the annealed 
mean- field predictions, must thus be explained in terms 
of the non-local constrain of absence of loops imposed 
in this kind of graphs, which is hard to implement in 
theoretical approaches. 

In this paper we explore the peculiarities induced in dy- 
namical processes by the absence of loops by considering 
the simplest possible example, namely the uncorrelated 
random walk [1?], OH • Several works have been devoted 
in the past to the study of random walks on complex 
networks, showing in general a good agreement between 
theory and simulations on looped networks, while differ- 
ences were reported in tree networks in Ref. [14J. Here, 
we find that the global constraint of lack of loops induces 
a general slowing down of diffusion, as measured by the 
network coverage and the mean topological displacement 



from the origin. As well, it profoundly alters the degree 
dependence of the mean-first passage time. This is due 
to the fact that the source-target distance is dominating 
in trees. In order to account for this features, we study 
the mean round trip time versus distance and find an 
analytic expression of its dependence on degree. 



II. RANDOM WALKS ON COMPLEX 
NETWORKS AND TREES 

We consider random walks on general networks defined 
by a walker that, located on a given vertex of degree 
k at time t, hops with probability 1/fc to one of the k 
neighbors of that vertex at time t + 1. We have mea- 
sured the properties of random walks on growing SF trees 
created with the linear preferential attachment (LPA) 
algorithm [ill [l9j : at each time step s, a new vertex 
with m edges is added to the network and connected 
to an existing vertex s' of degree k 8 i with probability 
lists' = (k s i + a) /(2m + a)s. This process is iterated 
until reaching the desired size N. The resulting network 
has degree distribution P(k) ~ fc -7 with tunable expo- 
nent 7 = 3 + a/m, with 7 < 3 for a < 0. For m = 1 
the LPA model yields a strict tree topology. Degree cor- 
relations, measured by the average de gree of the nearest 
neighbors of the vertices of degree k [20j, are given by 
knn(k) ~ 7V(3-7)/(7-i)fc-3+7 ^ . Therefore, only for 
7 = 3 (a = 0) we expect to obtain uncorrelated networks. 

In order to explore the intrinsic properties of a tree 
topology, disregarding SF effects, we have also consid- 
ered homogeneous networks. In the growing exponential 
network model (EM) [9J , at each time step s a new vertex 
with m edges is added to the network, and it is connected 
to m randomly chosen other vertices. In the continuous 
degree approximation (i.e., considering the degree as a 
continuous variable and substituting sums by integrals), 
this models leads to networks with an exponential degree 



distribution, P(k) — 



X—h/m 



jm. Again, homogeneous 



trees are generated by selecting m = 1. The random 
Cayley tree (RC), on the other hand, is generated by 
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adding z neighbors to a randomly selected leaf (i.e. a 
vertex whose degree is k = 1) at each time step s (z + 1 
neighbors are added to the first vertex). The resulting 
tree contains only vertices with degree k — z + 1, and 
leaves with k = 1. 

To check our results against looped structures, we 
have considered the uncorrelated configuration model 
(UCM) [13] , yielding uncorrelated networks with any pre- 
scribed SF degree distribution. The model is defined as 
follows: (1) Assign to each vertex i in a set of N ini- 
tially disconnected vertices a degree ki, extracted from 
the probability distribution P(k) ~ /c -7 , and subject to 
the constraints m < ki < A 1 / 2 and ^ ki even. (2) Con- 
struct the network by randomly connecting the vertices 
with edges, respecting the preassigned degrees 

and avoiding multiple and self-connections. Using this 
algorithm, it is possible to create SF networks which are 
completely uncorrelated. Additionally, by selecting the 
minimum degree m > 2, we generate connected networks 
with probability almost 1. The effect of correlations in 
looped structures can be checked by means of the config- 
uration model (CM) [§], which is analogous to the UCM, 
but allows degrees to range in the interval m < ki < N 
[23l ] . In all present simulations, we set for looped neworks 
to = 4, tree networks corresponding to m = 1 (z = 4 for 
the RC tree). 



III. RANDOM WALK EXPLORATION 

We start by studying two properties of a random walk 
that quantify the speed at which it explores its neigh- 
borhood in the network. The first one is the coverage 
S(t), defined as the number of different vertices visited 
by a walker at time t, averaged for different random walks 
starting from different sources. For looped networks, the 
coverage reaches after a short transient the functional 
form [24], [25| Sl (t) ~ t [43] , in accordance with theoreti- 
cal calculations for the Bethe lattice [26| , and eventually 
saturates to Sl (00) = N, due to finite size effects. A 
scaling form for the coverage has been proposed [24J to 
be S L (t) = Nf(t/N), with f(x) ~ x for x < 1 and 
f(x) ~ 1 for x > 1. 

The origin of the scaling of the coverage with system 
size can be understood by means of a simple dynamic 
mean-field argument. Let us define pk{t) as the proba- 
bility that a vertex of degree k hosts the random walker 
at time t. During the evolution of the random walk, this 
probability satisfies, in a general network with a correla- 
tion pattern given by the conditional probability P(k'\k) 
that a vertex of degree k is connected to another vertex 
of degree k' , the mean-field equation 



= -Pk(t) +k2_^ y Pk'{t). 



dt 



(1) 



In the steady state, d t pk{t) — 0, the solution of this equa- 
tion, for any correlation pattern, is given by the normal- 



10' 

io' J 
10- 4 

Iff' 

m> 
it)" 2 
'iff 4 
iff* 







CM y=2.5, k=5 ! 




— N=10 4 






"N=10 5 






— N=10 6 


1 













N=10 6 j 




. . CM y=2.5, k=5* , * b! * a -_ 






m — ■ CM Y=2.5, k=15 '- i 






— . UCM y=3.0. k=5 






~ UCM 1=3.0: k=15 






~ EM, k=5 






— EM, k=15 

1 





5 10 
kt / (<k>N) 



FIG. 1: (Color online) Coverage spectrum Sk(t) in looped 
complex networks as a function of kt/{k)N. Top: Curves 
for the same degree and different network size N. Bottom: 
Curves for different degrees and fixed network size. 
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Let us now define the coverage spectrum Sfc(i) as frac- 
tion of vertices of degree k visited by the random 
walker at least once. Obviously, we have that S(t) = 
NJ2k P(k)sk(t). The spectrum Sk{t) increases in time 
as the random walk arrives to vertices that have never 
been visited. Therefore, at a mean-field level, it fulfills 
the rate equation 



ds k (t) 
dt 



(3) 



Approximating pk> (t) by its steady-state value (for not 
too small times), we obtain 



ds k (t) . , k 

--[I- s k (t)\- 



dt l ^ " KX " n (k)N : 
whose solution, with the initial condition Sfc(0) = is 

Sfc(t) = 1 - exp ("(^v) 

We therefore are lead to the general scaling expression 
S(t) = ^ _ 

k 



(4) 
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In the limit of kt/(k)N 1, we recover the exact result 
S(t) ~ t (2(|. For SF networks, we obtain within the 
continuous degree approximation 



N 



7-1 



= 1 — (7 — 1)to 



exp 



kt \ 



dk 



(7) 



3 




10" 
t/N 



CMy=2.5. N=10 
CM y=2.5, N=10 5 
CM y=2.5, N=10 6 
UCM 7=3.0, N=10" 
UCM 7=3.0, N=10 3 
UCM 7=3.0, N=10 6 
EM,N=10 J 
EM, N=10 5 
EM, N=10 6 



10' 

io~ 2 
io- 4 

10° 
10" 
10" 
10" 





1 1 1 1 1 

LPAy=2.5,k=5 


' ' N=10 4 ^"^^ 




! — " N=10 3 




1 — N=10 6 




i 



. LPA 7=2.5. k=5 
■ LPA 7=2.5, k=15 
.LPA 7=3.0, k=5 
. LPA 7=3.0, k= 15 
EM Tree, k=5 
EM Tree, k=15 




300 



1011 



200 



300 
kt / (<k>N) 



400 



500 



dOO 



FIG. 2: (Color online) Rescaled coverage for looped com- 
plex networks. We plot in full lines the analytical predictions 
Eqs. ([7]) and JH), corresponding to SF and EM networks. Re- 
sults for CM, UCM, and EM networks have been shifted (top 
to bottom) in the vertical axis for clarity. 



FIG. 3: (Color online) Coverage spectrum Sk(t) in tree net- 
works as a function of kt/{k)N. Top: Curves for the same 
degree and different network size TV. Bottom: Curves for 
different degrees and fixed network size. The scaling here is 
different from the mean-field prediction Eq. ([5}. 



where E 1 {z) is the exponential integral function [2a ]. For 
EM networks, on the other hand, we find 
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In Fig. [TJ we can observe that the scaling predicted 
by Eq. ([5]) for the coverage spectrum Sk(t) is very well 
satisfied in looped complex networks, independently of 
their homogeneous or SF nature, and in this last case, 
of the degree exponent and the presence or absence of 
correlations. In Fig. [21 on the other hand, we plot the 
total coverage S(t) /N, which can be fitted quite correctly 
by the analytical expressions Eqs. ([7]) and ([H]) for SF and 
EM networks. 

On tree networks we find a different scenario. In Fig. [3] 
we can see that the coverage spectrum does not scale as 
predicted by our mean-field argument. While we do not 
have theoretical predictions for the correct scaling form, 
a numerical analysis of the total coverage, Fig. [U shows 
that, at short times, it grows in trees as Sx(t) ~ t/ln(t), 
preserving an approximate scaling form 



St® = N f 



t 



ln(t)N 



(9) 



with a scaling function f(x) that depends slightly on 
the network details (degree exponent, correlations, etc.). 
This observation indicates the presence of a general slow- 
ing down mechanism in the random walk dynamics in 
trees: the dynamics turns out to be more recurrent and 
therefore it is more costly to find new vertices during the 
walk. It is easy to see that this situation will correspond 
to a walker deep in the leaves of a subtree that has other- 
wise completely explored. In order to find new vertices, 
the walker must first find the exit to the subtree. This 
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FIG. 4: (Color online) Rescaled coverage as a function of time 
in complex trees. Results for LPA (7 = 2.5), LPA (7 = 3.0), 
EM, and RC networks have been shifted (top to bottom) in 
the vertical axis for clarity. 



difficulty in finding new vertices can be directly measured 
by the time lag At between the discovery of two new ver- 
tices. In Fig. \5\ we plot the probability distribution of 
time lags, P(At), computed for the discovery of the first 
1% of the network, for looped and tree structures. We 
observe that, in looped networks, this distribution takes 
an exponential form, compatible with an almost constant 
time lag between the discovery of two new vertices. In 
tree networks, this distribution shows instead long tails, 
that can be fitted to a lognormal form, indicating that, 
in some events (i.e. when the walker is trapped in one 
leaf in a subtree) the discovery of a new vertex can take 
an unusually large time. 

Acute signatures of slowing down can be found also in 
the analysis of the mean topological displacement (MTD) 
d(t) of the walker from its origin at time t, defined as the 
shortest path length from the vertex of origin of the ran- 
dom walk, to the vertex it occupies at time i, averaged 
over different source vertices, and different random walk 
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FIG. 5: (Color online) Distribution of lag times in tree (top) 
and looped (bottom) networks can be fitted, respectively, by 
a log-normal (full lines) and an exponential (dashed line) dis- 
tribution. Data refer to graphs of size N = 10 5 . 
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FIG. 6: (Color online) MTD as a function of time for looped 
complex networks of size N = 10 6 . 



realizations. In other works, the mean square topolog- 
ical displacement (RMSTD) d 2 (t) [H, [23] was instead 
considered. In complex networks, and since the shortest 
path length is a positive definite quantity, both quantities 
yield the same scaling result, i.e. d 2 {t) ~ [rf(t)] 2 [44j . 

In looped networks, see Fig.[6l we observe a very rapid 
increase of the MTD with time. In previous works [29(, 
the growth of the average distance (there measured in- 
stead as the RMSTD) was found to be a power law, 
d(t) ~ t a at early times in SF networks, with the ex- 
ponent a depending on the degree exponent. In our sim- 
ulations on looped networks generated with different al- 
gorithms, we do not find a clear signature for a power 
law behavior, which can only be approximately found 
in a tiny range of values of t at the very beginning of 
the walk, even for networks of size N — 10 6 . At large 
times, on the other hand, the MTD reaches a plateau, 
d(oo) — (d), due to finite size effects. The value of this 
plateau can be estimated using simple quantitative ar- 
guments. Assume that the random walker starts from a 
source vertex of degree k. During its dynamics in the 
steady state, it visits vertices of degree k' with probabil- 
ity, see Eq. @, k'/(k)N. On the other hand, vertices of 



FIG. 7: (Color online) Rescaled MTD as a function of time 
for looped complex networks. 



degree k and k' are, in average, at a topological distance 
[30l.l3ll] dk,k' ; therefore, we will expect the random walker 
to be at an average distance of a source of degree k 



(10) 



A further average over all possibles sources, leads to an 
average distance of the walker, for any source vertex, 
given by 



kP{k) 



(11) 



where dk = J2k> P(k')dk,k' is the mean topological dis- 
tance from any vertex to a given vertex of degree k [3ll | . 
The scaling of (d) with system size can be easily predicted 
assuming the expressions of dk in Ref. [3l[, namely 



A In 



N 



fc(7"l)/2 



(SF networks) 



(12) 



dk — A'lnN — B'k (exponential networks), (13) 

where A, A' and B' are size-independent constants. This 
yields in both cases 



(d) ~ InN. 



(14) 



Turning to the numerical data for looped networks in 
Fig. [3 we observe that it is compatible with a scaling 
behavior of the form 



-/,(/. = u/)/(-t 



(15) 



This scaling indicates that, after a short characteristic 
time t c ~ (d) ~ In N, the walker is in average as far as 
the origin at it can be, and it can therefore freely explore 
the whole network. 

In trees, Fig. [8j on the other hand, we observe a much 
slower growth of the MTD at early times, which can be 
approximately fitted with the form 



d T (t) ~ (In*)", 



(16) 
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FIG. 8: (Color online) MTD as a function of time for complex 
trees of size N — 10 6 . The exponent a is determined by a 
numerical fit from the observed relation dr(t) ~ (hit)" (data 
not shown). The values used in Figure are a ~ 0.91 for LPA 
trees with 7 = 2.5, a ~ 1.04 for LPA trees with 7 = 3.0, 
a ~ 1.31 for EM trees and a ~ 1.62 for the RC tree. 



where the exponent a depends on the details of the net- 
work. The whole function dx (t) is also observed to fulfill 
the scaling form 



d T (t) = (d)f 



hit 

(dy/c 



(17) 



This form implies that the characteristic time to es- 
cape from the neighborhood of the origin scales as t c ~ 
exp((d) 1 / Q ) ~ exp[(ln N) 1 ^], which means that the ex- 
ploration process is much more slower in trees, with the 
walker spending large amounts of time exploring the close 
vicinity of the origin of the walk. We remark that here 
the scaling function f{x) displays some further depen- 
dences on degree exponent, average degree and degree 
correlations in both looped and tree networks. 

The fact that the presence of a tree- like structure slows 
down the distance explored by a random walker on a net- 
work, allows to interpret the results presented in Ref . [2!| , 
in particular the power-law behavior at initial times of 
d(t). In fact, in Ref. 29] the substrate for the ran- 
dom walk simulations were SF networks generated with 
the CM model with minimum degree m = 1. In this 
case, simulations were performed on the giant compo- 
nent. Apart from the possible effect of degree correla- 
tions for 7 < 3, the point is that, for m = 1, traces 
of tree-like structure are still present in the network, in 
the form of chains of small degree vertices [Hj]. Thus, 
a remnant slowing down effect of the tree component is 
observed, see Fig. [SJ leading to an MTD that, at short 
times, scales as d(t) ~ t - 55 for the data at m = 1 shown 
in this graph, in excellent agreement with the observation 
in HI, namely d 2 (t) ~ t 11 for the RMSTD. 

A further remark concerns the relation between our 
results and the above mentioned analytical calculations 
for Bethe lattices [26| . according to which these struc- 
tures exhibit a behavior analogous to the one observed 
in looped networks. The apparent incongruity vanishes 




FIG. 9: (Color online) MTD as a function of time for UCM 
looped networks (7 = 3.0 and N = 10 6 ) with varying mini- 
mum degree m. 



when noticing that, while Bethe lattices are infinite hier- 
archical structures, we have focused on complex (i.e. dis- 
ordered) finite trees. To recover numerically the Bethe 
lattice behavior, indeed, it is necessary to adopt special 
algorithms in order to simulate an infinite hierarchical 
tree [3! HI. 



IV. MEAN FIRST-PASSAGE TIME 

More information about the dynamics of random walks 
can be extracted from the analysis of the mean first pas- 
sage time (MFPT) [35] r(i — > j), defined as the average 
time that a random walker takes to arrive for the first 
time at vertex j, starting from vertex i (27| . In networks 
with no translation symmetry, the MFPT from a source 
i to a target j needs not be equal to the MFTP from 
source j to target i. Therefore, different reduced MF- 
PTs can be considered. We can thus define the direct 
MFTP T~*(k) as the MFPT on a target vertex of degree 
k, starting from a randomly chosen source vertex, and 
the inverse MFPT r^(fc) as the MFTP on a randomly 
target vertex, starting from a source vertex of degree fc, 
namely 



r-(fc) 



jev(k) 



^( fc ) = ^E^ E 

j iev(fe) 



(18) 
(19) 



where V(k) is the set of vertices of degree k and Nk is 
the number of such vertices. A simple argument can 
predict the form of the MFPTs for random uncorrelated 
networks. In this case, the probability for the walker to 
arrive at a vertex i, in a hop following a randomly chosen 
edge, is given by q(i) = q(ki) = ki/(k)N (36[. Therefore, 
the probability of arriving at vertex i for the first time 
after t hops is P a (i;t) = [1 — q^)] 1 ^ 1 q{i) . The direct 
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FIG. 10: (Color online) Reduced MFPTs as a function of 
the degree k for looped complex networks. We recover the 
simple mean-field predictions r£*(k) ~ (k)N/k (main figure) 
and rr(fc) ~ (k){k- 1 }N (inset). 

MFTP to vertex i can thus be estimated as the average 

T^(ki) = ^2tP a (i',t) = ^-. (20) 

t 1 

For the inverse MFPT, we notice that, in a random net- 
work, after the first hop, the walker loses completely the 
memory of its source degree, therefore we can approxi- 
mate 

r^(fc) = £P(fc)r^(fc) = (fcXfc- 1 )^. (21) 

k 

Less trivial approaches [I?], H3] show in fact that the 
MFPT from a source vertex i to target vertex j depends 
on the degree of the target vertex as r(i — > j) ~ 1/kj, 
but has a residual dependence on the source vertex and 
it is actually asymmetric, r(i — > j) ^ r(j — ► i). This fact 
could in principle affect the form of the reduced MFPTs 
in real networks, defined in Eqs. (fl~g)l and (fT^]) . Fig. [TU1 
however, shows that for looped networks the behavior 
predicted for random uncorrelated networks turns out to 
be extremely robust with respect to changes in the topo- 
logical properties of the network: homogeneous or het- 
erogenous nature, degree exponent, presence or absence 
of correlations, etc. [23, [53, [III . 

In trees, on the other hand, we find a completely dif- 
ferent picture, see Fig. [TTJ Now, the direct MFPT in SF 
trees decays with k much slower than in looped networks. 
In fact, we can fit it numerically to the form 

T ^(k) = CiNlnN -C 2 Nln(k + C 3 ), (22) 

where C\ , Ci and C3 fitting parameters that depend only 
slightly on the network size. The NlnN dependence 
can be directly observed by plotting r^T*(l) for different 
system sizes, as shown in Fig. [T2J For homogeneous EM 
networks, on the other hand, the direct MFPT can be 
fitted to the form 

(fc) = D ± N InN- D 2 Nk, (23) 
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FIG. 11: (Color online) Direct MFPT as a function of the 
degree k for SF tree networks (iV = 10 4 ). Dashed lines cor- 
respond to nonlinear fittings to the empirical form Eq. ((22}. 
Inset: Direct MFPT as a function of the degree k for homo- 
geneous EM tree networks. The dashed line corresponds to a 
fitting to the empirical form Eq. (|23[1 . 

see inset in Fig. [TT] The scaling of t^*(1) in this case 
is also checked in Fig. [T3J With respect to the inverse 
MFTP, it is again constant, but now scales with system 
size as r^r(k) ~ NlnN for all kinds of trees (inset in 

Fig. HD. 

The topological structure of the trees can explain the 
unusual form of the MFPTs. While in looped networks 
the number of access paths to the target vertex is re- 
lated to its degree, on the tree the path is unique, and 
is given by the one-dimensional set of links and vertices 
connecting the starting vertex to the target. In this case, 
the degree of the target is much less important from the 
point of view of the walker, since finding the target cor- 
responds to finding a particular leaf (i.e. a k = 1 vertex) 
of the sub-tree the random walker is exploring. This 
observation suggests that while in looped networks the 
MFPT into a vertex is dominated by its degree (because 
the latter is related the multiplicity of the entry paths to 
the vertex), in trees the distance between the source and 
the target can be much more relevant, and thus induce a 
larger MFPT. 

We therefore consider the MFPT as a function of the 
topological distance dij between the starting vertex i and 
the target j [13, H^]. Since the distance between two 
vertices is by definition a symmetric quantity, it seems 
natural to re-define the MFPT in terms of the symmetric 
mean round trip time (MRTT) 

f(dy)=T(i-i) + TCj-»»). (24) 

i.e. the average time to go from i to 7 and back or vice- 
versa. It has been recently proved [33] that, for complex 
scale-invariant networks, the MRTT averaged for all ver- 
tices at the same distance scales as 

f(d) ~Nd D ™- D \ (25) 

where is the box dimension of the network, and D w 
its walk exponent [40( . For a class of scale- invariant net- 
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FIG. 12: (Color online) Direct MFPT on leaves in complex 
trees as a function of the network size TV. The observed scaling 
is 7"t*(1) ~ TV In TV. Inset: Inverse MFPT on complex trees 
for different network sizes (TV = 10 3 full colored points, TV = 
3 x 10 3 light colored points, TV = 10 4 empty points). The 
observed scaling is again rf"(fe) ~ TV In TV. 
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FIG. 13: (Color online) MRTT f(d) as a function of the 
source-target topological distance d in trees. Different curves 
collapse perfectly on fr(d) ~ Nd. 



works 40] corresponding to a tree structure, for which 
D w — Db = 1, the authors of Ref. [3!| obtained corre- 
spondingly a linear scaling fx (d) ~ Nd. We have checked 
that this linear form holds for different SF, EM and RC 
trees, see Fig. [T31 a result that leads us to conjecture 
that, for any complex tree, D w — = 1. 

We can use the result in Eq. (|25p to gain insight on 
the behavior of the anomalous reduced MFPTs in tree 
networks. Considering an average over all vertices with 
the same degree, we have that f(dkk') = ~~ * + 
r(k' — > fc). Averaging now over fc', we can consider the 
reduced MRTT 

f (fc) = ^ P(A:')[T-(fc ^ fc')+^(^' -» *)] = T^(k)+T^(k), 

k' 

defined as the average time to go from a randomly chosen 
vertex to a given vertex of degree k, and back (or vice- 
versa, since the MRTT is symmetric). Now, since f(dkk') 




FIG. 14: (Color online) Rescaled MRTT f T (t) as a function 
of the source degree k and for randomly chosen targets in SF 
trees. Predictions of Eq. (30), i.e. A^M/Mi)-! _ ^-t)/ 2 , 
are plotted as dashed lines. 



is linear in dkw for tree networks, we have 

f(fc) ~ ^ P(k')Nd kkl = Nd k . (27) 
fe' 

Assuming the scaling of d k as given by Eqs. (|T2)) and (fTB|) . 
we obtain 



f T (fe) ~ TV A In 



TV 



fc(7"l)/2 



for SF networks and 



(28) 



(29) 



f T {k) ~ TV (A' In TV - 73'fc) 



for EM networks. The unknown constant in Eq. 
can be reabsorbed in the the value of fr(l), to obtain a 
scaling form with system size for SF networks that reads 



vr(k) 



1 



In 



f T (l) In TV V> (7 " 1)/2 



TV 



(30) 



In Fig. [2] we show that this scaling form is very well 
satisfied by the MRTT in SF trees, independently of the 
degree exponent and correlation patterns, at least for in- 
termediate values of k. The observed bending at small 
degrees can be ascribed to the presence of a constant in 
the logarithm analogous to empirical parameter C3 in 
Eq. (f22]) . that does not follow from our argument. Finite 
size effects, on the other hand, are responsible for the 
deviations present at large degrees, that are indeed more 
evident in SF trees with smaller values of 7. 

This observations allow us to interpret the anomalous 
functional form of the reduced MFTPs observed in trees. 
From Eq. (|2l)|) , we have 

TT(k)=Tr(k)-Tr(k). (31) 

Writing r^(fc) ~ CTVlnTV, from Eq. (2SJ we obtain, for 
SF networks, 

(fc) ~ (A - C)N lnTV - A<yl ~ ^ TV In k (32) 
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while for homogeneous EM networks, we have 

(k) ~ (A' - C)N In N - B'Nk, (33) 

in agreement with the empirical fitting found in Eqs. (|22j) 
and pg]). 

This argument cannot be extended to looped networks, 
since here f{dkk') is not linear in dkk' ■ The k depen- 
dence of the MRTT can be however trivially obtained 
from the reduced MFPTs as f L (k) = T£*(k) + r^(k) ~ 
(k)N((k- x ) + l/k). 

V. CONCLUSIONS 

In this paper we have shown that complex tree-like 
topologies heavily affect the behavior of a random walk 
performed on top of them, with a global slowing down of 
the dynamics and a logarithmic dependence of the first 
passage time properties in SF networks. These features 
are intrinsically connected with the complex tree struc- 
ture and cannot be attributed to the mere presence of 
leaves, while they are radically different from the ones 
exhibited by Bethe lattices, i.e. infinite and hierarchical 
tree structures. 

We have studied the random walk exploration prop- 
erties and we have shown that complex trees induce 
a slower dynamics, compared to looped networks, for 



both the coverage and the mean topological displace- 
ment problems. Moreover, by means of the analysis of 
the symmetrized MFPT (the MRTT), we have been able 
to recognize the different role played by the degree k of 
the target vertex in looped and tree structures. In the 
former, a larger degree corresponds to a larger number 
of access ways to the target vertex. In the latter, on the 
other hand, the target vertex is always seen as a leaf by 
the random walker, and its degree k affects the MFPT 
only through the dependence of the average distance d^ 
between it and the rest of the vertices. These results pro- 
vide important insights into diffusion problems on trees, 
and help explaining the characteristic slow dynamics ob- 
served on diffusive processes taking place on top of tree 
networks [HI, EH, E2] • Moreover, they are also interesting 
in the study of dynamics in real- world networks, in which 
the so-called border trees motifs [4l[ have been recently 
shown to be significantly present . 
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